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Improving  the  efficiency  of  thermoelectric  devices  is  critical  to  their  widespread  adoption.  Here  a  design 
methodology,  formulated  on  computational  and  analytical  modeling,  derives  the  optimum  efficiency  and 
geometry  of  segmented  Bi2Te3-PbTe  Thermoelectric  Generators  (TEGs)  between  =298  K  and  =623  K 
(AT  =  325  K).  Comparisons  between  the  different  TEG  designs,  in  terms  of  the  electrical  load  to  TEG  elec¬ 
trical  resistance  ratio  (m  =  RL /RT£G),  are  simplified  thanks  to  the  devised  maximum  efficiency  to  temper¬ 
ature  gradient  (pmax  =  r;/AT)  metric.  Quasi-computational  results  of  [lmax  show  that  the  collective  Seebeck 
coefficient  Bi2Te3-PbTe  (a)  design  sustains  a  higher  electrical  load  in  relation  to  the  homogeneous  Bi2Te3 
and  PbTe  materials.  The  average  (a)  and  collective  (a)  Seebeck  coefficient  Bi2Te3-PbTe  configurations,  in 
comparison  to  Bi2Te3  and  PbTe,  exhibit  a  considerably  higher  (60-68%)  source  and  sink  thermal  resis¬ 
tance  matching  ( 0teg=  ®h*)-  The  proposed  segmented  Bi2Te3-PbTe  (a)  TEG  yields  a  peak  efficiency  of 
5.29%  for  a  AT  of  324.6  K. 

©  2012  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Photovoltaics  and  concentrated  solar  power  currently  dominate 
the  renewable  solar  power  generation  landscape.  In  meeting  the 
prevailing  and  imminent  energy  demands  of  the  world,  given  the 
comparative  advantages  of  each  technology,  it  will  be  critical  to 
diversify  the  power  generation  mix  with  the  integration  in  the 
building  and  transportation  sectors  of  other  solar  energy  conver¬ 
sion  technologies.  Enhanced  efficiency  and  reduced  cost  thermo¬ 
electrics  (TE),  which  directly  convert  thermal  energy  into 
electricity,  are  poised  to  make  significant  contributions  towards 
this  direction. 

Thermoelectric  (or  Seebeck)  devices  are  solid  state  electric  gen¬ 
erators  with  no  moving  parts.  Semiconductor  materials,  with  TE 
properties,  when  subject  to  a  temperature  gradient  across  their 
ends  produce  an  electromotive  force  (emf).  Efficient  TE  materials 
are  characterized  by  a  high  thermoelectric  figure  of  merit  (=1 ): 
z-  T,  where  z  =  ((a2«r)/ic),  T  is  the  absolute  temperature  (K),  a  the 
Seebeck  coefficient,  a  the  electrical  conductivity,  and  k  the  thermal 
conductivity.  Advantages  such  as  the  high  degree  of  reliability, 
simplicity,  tacit  operation,  lack  of  vibrations,  and  scalability  render 
Thermoelectric  Generators  (TEGs)  ideal  for  small  scale,  distributed 
power  generation,  and  energy  harvesting  applications  [1]. 
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According  to  the  compatibility  factor,  herein  denoted  by 

i-(Vi+z-r-i)/{«r)  (i) 

the  efficiency  of  segmented  TEG  (STEG)  materials  increases  when  the 
two  joined  TE  materials  are  compatible,  that  is,  their  s  values  do  not 
differ  by  a  factor  of  two  or  more  [2],  Application  of  the  compatibility 
factor  to  segmented  and  cascaded  TEG  designs  showed  that  the  use  of 
incompatible  TE  materials  results  in  a  drop  in  efficiency  [3], 

Progress  in  sophisticated  nanoscale  TE  materials  as  well  as 
strategies  to  enhance  their  thermopower  and  reduce  their  thermal 
conductivity  have  been  outlined  by  Snyder  et  al.  [4],  Other  investi¬ 
gators,  using  computational  tools  and  a  prototype  model,  have 
examined  the  influence  of  the  heat  exchangers’  thermal  resistance 
on  the  TEG  efficiency  [5],  Depending  on  the  particular  circum¬ 
stances,  the  application  of  TEGs  on  vehicles  for  heat  recovery  re¬ 
sulted  in  unfavorable  and  beneficial  results  [6],  More  recently,  a 
flat-panel  solar  thermoelectric  generator  device,  incorporating 
high-performance  nanostructured  TE  materials  and  spectrally- 
selective  solar  absorbers,  in  an  evacuated  environment,  with  a 
demonstrated  maximum  efficiency  of  4.6%  has  been  proposed  [7], 
One  of  the  most  effective  ways  to  boost  the  efficiency  of  a  thermo¬ 
electric  generator  is  to  establish  a  large  temperature  gradient  (AT) 
across  the  TEG  constituent  elements.  In  applications  other  than  ther¬ 
mal  energy  waste  recovery,  concentrated  solar  energy,  using  rela¬ 
tively  inexpensive  optics  (e.g.,  Fresnel  lens),  can  help  create  a  large 
temperature  differential  across  a  TEG.  Noticeably,  even  with  a  high 
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Nomenclature 

Ay 

cross-sectional  area  of  p-  and  n-type  TEG 

Z 

device  thermoelectric  figure  of  merit  (K  ') 

thermoelements  (m2) 

E 

electric  field  (NT1) 

Greek  letters 

h0 

total  enthalpy  (m  2  s  1 ) 

a 

thermoelectric  material  Seebeck  coefficient  (V  1<  ’) 

I 

electric  current  (A) 

a 

average  TE  material  Seebeck  coefficient  (V  K  1 ) 

im  ox 

maximum  current  density  (A  m  2) 

a 

collective  TE  material  Seebeck  coefficient  (VK  ') 

i 

length  of  the  TEG  elements  (m) 

P  —  lr 

efficiency  to  temperature  gradient  ratio  (K  ’) 

m  =  ^c 

electrical  load  resistance  to  TEG  electrical  resistance  ra¬ 

n 

average  Bi2Te3-PbTe  Seebeck  coefficient  efficiency 

tio 

n 

collective  Bi2Te3-PbTe  Seebeck  coefficient  efficiency 

tl  =  Fhu 

TEG  to  heat  sink  and  heat  source  thermal  resistance  ra¬ 

0Hx 

heat  sink  and  heat  source  thermal  resistance  (KW1) 

tio 

&TEG 

TEG  thermal  resistance  (KW  ’) 

rlTEC 

number  of  thermoelement  pairs 

K 

thermal  conductivity  (W m-1  K  ’) 

P 

static  pressure  (N  nr2) 

A 

thermal  conductance  (WT1) 

Rl 

electrical  load  resistance  (C2) 

P 

dynamic  viscosity  (kg  m-1  s  ’) 

Rtec 

TEG  electrical  resistance  (C2) 

P 

electrical  resistivity  (C2  m) 

Tc 

heat  sink  temperature  (K) 

Pm 

mass-density  (kg  m  3) 

Th 

heat  source  temperature  (K) 

a 

electrical  conductivity  (£2  1  nr1) 

T Inter 

Bi2Te3-PbTe  interface  temperature  (K) 

r,j 

viscous  stress  in  the  j  and  t-directions,  respectively 

T inter 

average  Bi2Te3-PbTe  interface  temperature  (K) 

(N  nr2) 

Z 

material  thermoelectric  figure  of  merit  (K1) 

0 

electric  potential  0  C  ’) 

temperature  difference,  the  performance  of  a  TEG  varies  appreciably, 
largely,  due  to  a,  a,  and  k  being  temperature  dependent. 

In  circumventing  the  latter  limitation  two  or  more  TE  materials 
can  be  joined  in  series  to  form  a  segmented  TEG.  The  rationale  be¬ 
hind  the  STEG  design  is  to  boost  the  overall  TEG  system  efficiency 
by  having  individual  TEG  segments  operating  at  their  optimum  (or 
near-optimum)  z  ■  T  within  the  corresponding  operating  tempera¬ 
ture  range.  Augmented  TEG  efficiencies  have  been  recently  realized 
at  the  experimental  level  [8]  including  designs  for  space  applica¬ 
tions  [9], 

Here  a  design  methodology  is  proposed,  based  upon  computa¬ 
tional  and  analytical  modeling,  of  a  STEG  with  a  heat  sink  (cold) 
temperature  (Tc)  of  298.2  K  and  a  heat  source  (high)  temperature 
(Th)  of  622.8  K  (AT  =  324.6  K).  The  efficiency  of  bismuth  telluride 
(Bi2Te3),  lead  (II)  telluride  (PbTe),  and  the  average  (a)  and  collec¬ 
tive  (a)  Seebeck  coefficient  segmented  Bi2Te3-PbTe  configurations 
are  determined  as  a  function  of  the  electrical  load  resistance  to  the 
TEG  electrical  resistance  ( m  =  RJRtec )  ratio  and  the  TEG  and  heat 
source  and  sink  thermal  resistance  matching  (  0Tec  =  &h*)- 

The  segmented  Bi2Te3  (ot)  system  overestimates  the  maximum 
efficiency  to  temperature  gradient  ratio  (/imax  =  rj/AT)  while  Bi2_ 
Te3-PbTe  (a)  yields  a  [!max  value  which  resides  intermediate  to 
the  pmax  values  of  Bi2Te3  and  PbTe.  Technically,  it  is  possible  to  de¬ 
velop  the  preceding  temperature  gradient  of  AT  325  K  using 
inexpensive  Fresnel  lens,  at  a  concentration  factor  CF  »  50,  where 
CF  is  the  ratio  of  the  optically  active  collector  area  to  the  absorber 
area  upon  which  solar  radiation  impinges.  The  application  of  con¬ 
centrator  solar  optics,  which  serves  to  emphasize  such  a  possibil¬ 
ity,  has  been  extensively  covered  elsewhere  [10].  This  paper 
focuses  on  the  application  of  segmented  TE  materials  as  facilitated 
by  solar  thermal  energy. 


2.  Segmented  TEGs 

Considering  the  behavior  of  the  figure  of  merit  (z  ■  T)  for  TE 
materials,  as  a  function  of  the  operational  temperature  range  of 
622.8  K  and  298.2  K,  it  is  evident  that:  (a)  no  single  material  is 
thermally  stable  within  the  entire  temperature  regime  [11]  and 
(b)  higher  overall  performance  can  be  realized  if  two  (or  more) 
TE  materials  are  utilized  for  each  thermoelement  [12],  In  the  low 
temperature  range,  Bi2Te3  was  selected  due  to  its  high  figure  of 


merit  whereas  PbTe,  due  to  its  thermal  stability,  is  more  suitable 
for  the  top  temperature  regime.  Other  combinations  of  materials 
can  also  be  used,  if  desired.  Unless,  the  gain  in  system  performance 
is  justified  by  the  use  of  more  than  two  TE  materials,  connected  in 
series,  undue  structural  complexity  due  to  material  discontinuities 
-  resulting  also  in  higher  manufacturing  costs  -  is  something  to  be 
avoided. 

In  ascertaining  the  merits  and  shortcomings  of  various  design 
configurations,  different  theoretical  scenarios  were  investigated 


with  an  active  cooling  heat  exchanger  (figure  not  to  scale). 
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Table  1 

Lead  telluride  (PbTe)  TE  material  particulars. 


Parameter 


Peak  TEG  power  ( Wmp ) 
Max.  power  voltage  (Vrap) 
Seebeck  coefficient,  a 
PbTe  (n-type) 

PbTe  (p-type) 
Thermal  conductivity,  k 
PbTe  (n-type) 

PbTe  (p-type) 
Electrical  resistivity,  p 
PbTe  (n-type) 

PbTe  (p-type) 

Scale  sun  factor,  CF 


-225  pV  K  1 
+287  pV  K  1 

2.07Wm-,r' 
1.448  Wm'r1 

1.818  x  10-5£2m 
2.56  x  l(T5nm 


Reference 


Table  2 

Bismuth  telluride  (Bi2Te3)  TE  material  particulars. 


Parameter 


Peak  TEG  power  (Wmp) 
Max.  power  voltage  (Vmp) 
Seebeck  coefficient,  a 
Bi2Te3  (n-type) 

Bi2Te3  (p-type) 
Thermal  conductivity,  k 
Bi2Te3  (n-type) 

Bi2Te3  (p-type) 
Electrical  resistivity,  p 
Bi2Te3  (n-type) 

Bi2Te3  (p-type) 

Scale  sun  factor,  CF 


3  W 
5  V 

250  pVK  1 
+250  pVK  1 

1.359  Wm1  K- 
1.736  Wm-'T 

2  x  10  5  n  m 
2  x  10  5  £2  m 
1000  Wm2 


Reference(s) 


[15] 

[15] 

Calc,  from  [15,16] 
Calc,  from  [16] 

[16] 

[16] 


under  the  same  operating  conditions.  To  facilitate  our  analysis,  a 
design  methodology  was  developed  based  on  the  principles  of 
thermoelectricity  [11],  heat  transfer,  and  the  electrical  and  thermal 
characteristics  of  the  device,  the  electrical  load,  the  heat  sink  and 
the  heat  exchanger. 

As  depicted  in  Fig.  1,  the  segmented  thermoelectric  generator 
consists  of  a  pair  of  p-  and  n-type  Bi2Te3-PbTe  thermoelements 
joined  thermally  in  parallel  and  electrically  in  series.  Thermal  en¬ 
ergy  dissipation  of  the  TEGs  is  achieved  using  active  cooling  with 
water  as  the  cooling  medium.  Section  3.2. 1.2  outlines  the  heat  sink 
design  methodology.  The  Seebeck  coefficient  of  the  n-type  and  p- 
type  PbTe  TE  segments,  clustered  in  Table  1,  corresponded  to  dop¬ 
ing  levels  of  0.03  mol%  lead  (II)  iodide  (Pbl2)  and  0.3  atm.%  sodium 
(Na),  respectively.  In  the  case  of  n-  and  p-type  Bi2Te3  the  a  values, 
which  appear  in  Table  2,  refer  to  an  electrical  conductivity  (c)  of 
5  x  104  C2  1  nr1. 

Using  Fresnel  optics  it  is  possible  to  concentrate  photons  on  the 
top  plate  of  the  TEG  and  raise  the  heat  source  temperature  (TH)  to 
623  K.  Discrete  p-  and  n-type  TEG  legs  are  electrically  and  ther¬ 
mally  connected,  on  top  of  the  TEG,  using  a  copper  plate,  while  a 
heat  exchanger  on  the  lower  side  of  the  system  facilitates  thermal 
energy  dissipation.  In  contrast  to  natural  convective  heat  dissipa¬ 
tion,  forced  cooling,  as  explained  in  Section  3.2.1,  could  be  em¬ 
ployed  to  attain  more  effective  thermal  energy  discharge  at  the 
heat  sink.  System  particulars  are  presented  in  Tables  1  and  2. 

3.  Design  methodology  for  segmented  TEGs 

3.1.  Analytical  model 


the  latter  need,  in  this  work  the  formulation  of  a  practical  design 
methodology  for  segmented  TEGs  is  presented. 

Collectively,  the  TEG  is  to  produce  a  peak  power  ( Pmp )  of  8  W,  at 
a  maximum  power-voltage  ( Vmp)  of  1 1  V.  Because  the  TEG  could  be 
powered  by  concentrated  solar  thermal  energy,  system  efficiency 
is  investigated  between  a  peak  temperature  of  622.8  K  and  an 
ambient  temperature  of  298.2  K.  Maximum  allowable  current  den¬ 
sity  Umax )  was  set  to  1 00  A  cm  2  [11]  while  the  TEG  open-circuit 
voltage  ( Voc )  was  obtained  by  summing  the  Seebeck  coefficients 
(a)  of  the  constituent  p-  and  n-type  materials  over  the  temperature 
gradient  of  interest 

AVoc  =  j\<hPV)  -  *,,n(T)}dT  (2) 

where  i  is  the  material  of  interest  for  the  p-  and  n-type  materials, 
and  Th  and  Tc  are  the  source  and  sink  temperatures,  respectively. 
An  average  a  value  over  the  temperature  range  of  interest,  denoted 
by  a,  is  used  for  individual  Bi2Te3  and  PbTe  materials  and  the  seg¬ 
mented  Bi2Te3-PbTe  designs.  Given  that  the  TE  elements  are  con¬ 
nected  electrically  in  series,  the  (maximum  power)  current 
crossing  each  TEG  unit  can  be  determined  from  Pmp  =  Imp  x  Vmp. 
To  determine  the  number  of  TEG  units  (nrec)  needed  to  develop  a 
total  voltage  of  11  V,  the  design  methodology  proposed  by  Da  Rosa 
[11]  was  adapted  with  the  introduction  of  appropriate  modifica¬ 
tions  to  reflect  the  peculiarities  of  the  system.  The  maximum  power 
voltage  [ymp)  of  the  TEG  system  is  obtained  by  summing  the  open 
circuit  voltage  drop  across  resistance  R 

n  ■  Voc  -  n  ■  R  ■  I  =  Vmp  (3) 

To  deduce  the  number  of  TEG  pairs  (Hteg),  the  TEG  resistance  (R) 
needs  to  be  determined  first.  For  peak  efficiency,  the  load  resis¬ 
tance,  Rl,  is  given  by 

Rl  =  m  ■  Utec  ■  R  (4) 

where  m  =  y?i+  {T)  Z  and  (T)  =  (TH  +  Tc)/2.  Furthermore,  the  de¬ 
vice  figure  of  merit  (Z)  is 


where  A  and  R  are  the  thermal  conductance  and  the  electrical  resis¬ 
tance,  respectively.  The  product  A  ■  R  is  obtained  from 


where  i  and  j  denote  the  materials  of  interest  and  p  their  electrical 
resistivity.  Once  Z  is  found,  the  total  number  of  TEG  pairs  can  also 
be  obtained  at  the  corresponding  thermal  gradient  of  interest.  Con¬ 
ductive  heat  transfer  (PL)  through  the  TEG  is  given  by  PL  =  A  ■  ( TH  - 
-  Tc).  Under  electrical  loading  the  thermal  power  supplied  to  the 
TEG  by  the  heat  source  is 

PH  =  PL  +  aTHIi-^Rilf  (7) 

where  I  is  the  current  of  the  TE  material  of  interest  denoted  by  i. 
Considering  the  thermal  power  under  load  (PL)  and  no  load  (PH) 
conditions,  the  efficiency  of  single-material  Bi2Te3  and  PbTe  ther¬ 
moelements  can  be  determined  from 


(8) 


First  order  approximations,  based  on  simple  equations,  are  of¬ 
ten  used  to  arrive  to  approximate  design  characteristics  of  TEGs 
before  proceeding  into  more  rigorous  mathematical  methods,  such 
as  computational  modeling,  and  finally  experimental  fabrication. 
Hence,  reliable,  versatile  and  time  efficient  design  methodologies 
are  vital  to  the  progress  of  high  performance  TEGs.  Motivated  by 


The  average  Seebeck  TEG  design  efficiency  ( fj )  is  obtained  from 
the  mean  Seebeck  coefficient  (a),  in  turn,  derived  from  the  average 
total  Seebeck  values  of  the  Bi2Te3  and  PbTe  materials 
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Finally,  the  collective  TEG  design  efficiency  ( r j)  is  defined  by 
summing  the  individual  power  developed  by  each  of  the  constitu¬ 
ent  Bi2Te3  and  PbTe  thermoelectric  elements  based  upon  their  total 
Seebeck  coefficients 


(10) 


where  subscripts  i  and  j  are  the  thermal  power  generated  by  the  Bi2_ 
Te3  and  PbTe  thermoelectric  materials,  respectively. 

One  of  the  most  effective  and  yet  difficult  ways  to  improve  the 
efficiency  of  a  TEG  is  to  increase  the  device  figure  of  merit,  Z  (Eq. 
(5)).  Besides  selecting  TE  materials  with  the  largest  possible  z  val¬ 
ues,  the  geometry  of  the  TE  can  help  improve  performance  by 
decreasing  the  magnitude  of  the  denominator  (A  ■  R )  (Eq.  (6)).  It 
turns  out  that  there  is  a  specific  geometry,  in  terms  of  the  TEG 
length  and  cross-sectional  area,  which  minimizes  A  ■  R  [11] 


£i  Aj  /Kj  •  pA 1/2 

\KrPi) 


(11) 


where  t  is  the  length,  A  the  cross-sectional  area,  and  i  and  j  are  the 
p-  and  n-type  TEG  materials,  respectively.  Assuming  that  the  TEG 
arms  have  the  same  overall  length  and  experience  a  maximum  cur¬ 
rent  density  (Jmax)  of  100  A  cm  2,  A,-  and  Aj  can  be  determined.  Con¬ 
sidering  the  resistance  of  individual  thermoelements,  one  can 
obtain  the  length  of  the  TE  pairs  using 

R  =  Prli  +  Pi'Xj  w 

In  the  sequel,  the  performance  of  individual  Bi2Te3  and  PbTe 
materials  and  the  segmented  Bi2Te3-PbTe  designs  are  presented. 
Optimum  system  performance  was  obtained  in  the  context  of  the 
TEG  and  load  electrical  resistance  ratio  and  thermal  TEG  to  heat 
sink  and  exchanger  resistance  ratio.  Moreover,  the  optimum  geom¬ 
etry  for  each  TEG  configuration  was  determined.  Hereafter,  the 
computational  model  is  presented  and  the  integration  of  results 
between  the  analytical  and  the  computational  models  is  outlined. 


3.2.  Computational  Bi2Te3-PbTe  TEG  Model 

Computational  modeling,  through  the  solution  of  partial  differ¬ 
ential  equations  (PDEs),  generates  insight  otherwise  unfathomable 
by  experimental  means.  The  temperature  distribution  and  Direct 
Current  (DC)  conduction  across  the  thermoelements  which  also 
govern  the  TEG  efficiency  are  prominent  examples.  For  segmented 
TEGs,  of  particular  interest  is  the  temperature  gradient  across  each 
TE  material  which  is  dictated  by  the  material’s  boundary  condi¬ 
tions.  Knowledge  of  the  temperature  at  the  Bi2Te3-PbTe  interface 
can  be  used  as  input  in  the  analytical  model  to  derive  the  efficiency 
and  the  physical  dimensions  of  the  TEG  configurations  [17],  Joule 
heating  and  the  Peltier  and  Seebeck  effects  were  also  implemented 
in  the  computational  model.  Fluid  flow  in  the  heat  exchanger  con¬ 
duit,  as  shown  in  Fig.  2,  was  governed  by  the  solution  of  the  conti¬ 
nuity  and  the  Navier-Stokes  equations.  Table  3  presents  the 
physical  properties  of  the  solid  and  liquid  media  comprising  the 
computational  model. 

Fig.  2  shows  the  computational  mesh  used  to  trace  the  temper¬ 
ature  variations  across  each  TE  material.  The  computational  domain 
of  the  TEG  model  was  discretized  using  structured  conforming 
meshes  mainly  owing  to  the  relatively  simple  geometry  of  the 
TEG  design  and  the  higher  accuracy  permitted  by  a  structured,  ver¬ 
sus  an  unstructured,  grid.  The  computational  model  presented  in 
this  work  was  developed  using  the  multi-physics  CFD-ACE  +  plat¬ 
form  (ESI,  Paris,  France)  [18],  A  steady  state  time  dependence  was 
adopted  for  the  solution  of  the  governing  equations. 


6.0E-02 


Fig.  2.  The  TEG-heat  exchanger  computational  mesh  for  Bi2Te3-PbTe.  The  mesh 
consists  of  nine  distinct  computational  zones  with  a  total  of  39,069  structured  cells. 


Table  3 

Bi2Te3-PbTe  computational  model  particulars. 


Parameter  (units) 

Mass-density,  pm  (kg  m  3) 

Bi2Te3  (p-  and  n-type) 

PbTe  (p-  and  n-type) 

Copper (Cu) 

Stainless  steel  310 
Water  (H20) 

Specific  heat  capacity,  Cp  (J  kg~'  K  ’) 
Bi2Te3  (p-  and  n-type) 

PbTe  (p-  and  n-type) 

Stainless  steel  310 
Water 

Thermal  conductivity,  k  (Wm  1  K  ’) 
Copper 

Stainless  Steel  310 
Dynamic  viscosity,  p  (kg  m  1  s  ’) 


Value 


7740 


7782 

997 


186 

155 

385 

512 

4181 

401 


0.607 

8.94  x  10-4 


Reference(s) 


[20] 

[13] 


[211 


3.2 A.  Governing  equations  and  boundary  conditions 
3.2.1  A.  Heat  transfer.  Heat  flux  in  the  solid  domain  was  traced 
through  heat  conduction  while  thermal  energy  transfer,  in  the 
cooling  fluid,  was  governed  by  the  numerical  solution  of  the  two- 
dimensional  (2D)  total  enthalpy  equation  [22] 

d(uho)  d(vho)  =  (tfT  cZt\  (dim**)  8(UTyx)  d(tnw)  a(tn„)\ 

“TS~+Tf  +  v  at  +  #  +  ax  +  ay  ) 

(13) 

where  h0  =  i  +  p/pm  +  \  ■  (u2  +  v1)  is  the  total  enthalpy,  p  the  static 
pressure,  i  the  internal  energy,  and  tw  is  the  viscous  stress  with  i 
and  j  being  the  stress  components  acting  on  the  ^-direction  on  a  sur¬ 
face  normal  to  direction  i. 
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3.2. 1.2.  Heat  exchanger.  Establishing  a  high  temperature  gradient 
across  the  TEG  calls  for  the  use  of  an  effective  heat  sink  capable 
of  handling  a  large  thermal  flux.  For  this  purpose  we  have  investi¬ 
gated  both  passive  and  active  cooling.  Owing  to  the  high  thermal 
flux  a  heat  exchanger  needs  to  dissipate,  natural  convective  cooling 
-  due  to  the  large  length  of  the  cooling  fins  -  becomes  impractical. 

Energy  balance  calculations  revealed  that  heat  dissipation  can 
be  efficiently  accomplished  using  forced  cooling  with  water  as 
the  working  medium.  For  the  sake  of  simplicity  the  length  of  the 
conduit,  shown  in  Fig.  2,  was  assumed  to  match  the  same  axial 
span  as  the  thermal  bridge  connecting  the  two  TEG  arms.  Flow  in 
the  conduit,  in  Cartesian  coordinates,  is  governed  by  the  2D  conti¬ 
nuity  and  the  Navier-Stokes  equations  for  an  incompressible  New¬ 
tonian  fluid.  Considering  that  the  mass-density  ( pm )  of  the 
working  medium  does  not  vary  with  time  and  space,  the  2D  equa¬ 
tions  of  flow  with  constant  viscosity  [ft)  at  a  horizontal  pipe  orien¬ 
tation,  ignoring  gravity,  reduce  to  [23] 


du  dv_ 
dx  +  dy  ~ 

[x  direction) 
f  du 

P™(Udx+V 
(y  direction) 
f  dv 

Pm(unx  +  V 


du\ 

dy) 


dv\ 

'ay) 


dp  ( cPu  c^lA 

I+P{d*+W) 


_  dp  (d*V  d*V 

~~dy  +  p\dx2  +  dy2 


(14) 


(15a) 


(15b) 


3.2. 1.3.  Electric,  thermal  and  thermoelectric  effects.  Direct  current 
conduction  in  the  TE  model  was  implemented  through  the  numer¬ 
ical  solution  of  Laplace’s  equation  of  the  electric  potential  (<£).  The 
two-dimensional  Laplace’s  equation  is 


[  dx2  +  dx2  ) 


(16) 


where  a  is  the  electrical  conductivity.  Subsequently,  the  electric 
field  (E)  is  determined  from  the  2D  vector  field 


Considering  the  Seebeck  effect  (a),  Eq.  (16)  assumes  the  form 

(18) 


,»  (d2T  d2f 

*  -‘■’’[W+Wy, 


Similarly,  factoring  in  the  thermoelectric  phenomenon  the  new 
electric  field,  based  on  Eq.  (17),  is 


Joule  heating  emanating  from  the  thermal  energy  generated  in  a 
conductor,  due  to  current  flow,  was  modeled  as  a  heat  source  [18] 


Ht  =  a  ■  e*2  (20) 

Finally,  the  Peltier  power  (Ps)  implemented  at  the  interface  be¬ 
tween  two  materials  with  thermoelectric  properties  is  obtained 


Ps  =  a.  T  I  (21) 

where  /  is  the  electric  current. 


3.2.1. 4.  Boundary  conditions.  For  simplicity  the  top  copper  surface 
of  the  TEG  was  kept  isothermal  to  623  K  throughout  the  simula¬ 
tions.  The  rest  sides  of  the  model  were  set  adiabatic  permitting  thus 
only  heat  flux  through  the  TEG  elements  to  the  cooling  conduit  - 
the  latter  of  which  serves  as  the  heat  sink  (Fig.  2).  Water  entering 
the  conduit,  at  the  left  side,  was  assumed  to  possess  a  uniform 
velocity  of  0.1  m  s  1  at  a  constant  temperature  of  298  K.  Assuming 
a  pipe  with  smooth  internal  surface,  the  low  inlet  velocity  of  water 
generated  a  laminar  flow  (Reynolds  number,  Re  =  3,390).  Based  on 
this  laminar  flow  we  assumed  that  the  entrance  length  of  the  pipe 
does  exert  an  appreciable  effect  on  convective  heat  dissipation. 
Therefore,  we  did  not  consider  a  longer  length  pipe  section  neces¬ 
sary  for  establishing  a  fully  developed  flow. 

Current  density  at  the  top  electrical  bridge,  connecting  the  two 
TEG  elements,  was  set  to  100  A  cm2  while  the  bottom  ends  of  the 
TEG  legs  were  electrically  grounded  (</>  =  0).  No  current  was  as¬ 
sumed  to  escape  the  rest  sides  of  the  TEG  and  the  bottom  conduct¬ 
ing  bridge.  Hence  the  bottom  electrical  plate  and  the  cooling  pipe 
were  kept  in  electrical  isolation  from  the  rest  TEG  system. 


3.2.2.  Grid  dependence  study 

A  grid  (in)  dependence  study  serves  to  assess  the  sensitivity  of 
the  numerical  solution  to  parameters,  such  as,  the  type  and  number 
of  computational  cells.  Moreover,  it  offers  a  first  indication  as  to  the 
self-consistency  of  the  modeling  results.  Valuable  information  can 
also  be  extracted  regarding  the  trade-offs  between  the  accuracy 
and  reliability  of  results,  computational  expense,  and  time. 

To  raise  confidence  in  the  results  of  the  TEG  system,  the  temper¬ 
ature  distribution  of  three  meshes  consisting  of  19,852  cells  (Grid 
A),  39,069  cells  (Grid  B)  and  83,347  cells  (Grid  C)  were  compared. 
Temperature  evolution  was  traced  at  the  axial  center  of  the  n-  and 
p-type  rods  at  the  Bi2Te3-PbTe  junctions.  Converged  temperature 
results,  at  the  interface  of  the  n-  and  p-type  materials,  exhibited 
a  difference  of  less  than  1%  for  the  same  time  period. 

As  the  number  of  cells  increased  from  Grids  A  to  B  to  C  the  dif¬ 
ference  in  their  mean  interface  temperatures  (7W)  shrank.  Grid  B, 


T.r 


Stainless  Steel 


Fig.  3.  Temperature  distribution  in  the  Bi2Te3-PbTe  segmented  TEG,  at  t=  213  min. 
The  heat  source  (top)  temperature  (TH)  is  623  K  while  the  heat  sink  (bottom) 
temperature  (Tc)  is  298  K.  Minor  temperature  discrepancies  between  the  n-type 
(left)  and  the  p-type  (right)  TE  legs  arise  from  different  material  properties. 
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with  39,069,  was  selected  as  the  most  suitable  model  because  it 
exhibited  the  best  balance  in  terms  of  accuracy  and  computational 
time.  With  only  a  minute  T /ntEr  difference  of  3.5  x  10  2%  it  reached 
convergence  in  213  min  versus  351  min  for  Grid  C.  Coarser  grid  A 
was  rejected  because  it  yielded  a  TInter  deviation  of  0.88%  and 
0.92%  from  Grids  B  and  C,  respectively.  Section  3.2.3  outlines  the 
process  of  obtaining  TInter. 

3.2.3.  Computational  Bi2Te3-PbTe  TEG  model  results 

Numerical  solution  of  the  energy  equation  (13)  in  the  computa¬ 
tional  domain  (Fig.  2)  revealed  a  slight  temperature  discrepancy 
between  the  p-  and  n-type  interfaces  of  the  Bi2Te3  and  PbTe  TE 
materials,  as  depicted  in  Fig.  3.  To  avoid  undue  complexity  analyt¬ 
ical  calculations,  as  outlined  in  Section  4.2,  are  based  on  Tinter  ob¬ 
tained  from  the  computational  model.  Variations  in  temperature 
between  the  two  TEG  legs  emanate  from  differences  in  the  values 
of  the  physical  properties  of  the  thermoelectric  Bi2Te3  and  PbTe 
materials,  are  shown  in  Table  3. 


4.  Results  and  discussion 

4.1.  Load-TEG  electrical  resistance  matching 

To  evaluate  the  thermoelectric  compatibility  between  PbTe  and 
Bi2Te3,  the  compatibility  factor  (Eq.  (1))  was  used.  If  the  difference 
in  s  values  between  two  joined  materials  is  less  than  a  factor  of 
two,  the  collective  maximum  efficiency  does  not  deteriorate  with 
segmentation  [2],  A  temperature  gradient  of  153.2  K  across  Bi2Te3, 
obtained  from  the  computational  model,  yields  a  mean  s  value  of 
1.89  V  ’.  Lead  telluride,  within  the  temperature  range  of  operation 
(AT  =171.4  K),  generates  a  mean  s  of  1.602  V  '.  Comparing  the  lat¬ 
ter  two  s  values  reveals  very  good  thermoelectric  matching  for  Bi2_ 
Te3  and  PbTe.  Interconnected  Bi2Te3  and  PbTe  TE  elements 
including  the  copper  bridges  between  the  TE  legs  increase  the 
overall  electrical  resistance  of  the  TEG  system.  However,  for  small 
contact  resistances  (Rc)  in  relation  to  RTeg,  such  that  Rc  -c  Rtec.  re¬ 
sults  (not  included  here)  have  shown  that  the  contact  resistance 
does  not  appreciably  adversely  affect  the  efficiency  of  the  Bi2Te3, 
PbTe,  Bi2Te3-PbTe  (a),  Bi2Te3-PbTe  (a)  design  configurations. 

A  thermoelectric  generator  produces  maximum  power  output 
( Pmax )  when  the  TEG  electrical  resistance  (Rtec)  is  the  same  as  the 
load  electrical  resistance  (Rl)',  the  dimensionless  ratio  of  which  is, 
often,  denoted  by  m  =  RJRtec-  Eq.  (4)  offers  an  alternative  way  of 
deriving  m  but  not  in  terms  of  system  optimization.  Analyzing 
the  relationship  between  RL  and  i?TEG  helps  to  quantify  the  opti¬ 
mum  TEG  (power)  performance.  In  this  context,  there  exists  a  sin¬ 
gle  optimum  mopt  value,  usually  mopt  >  1,  at  which  the  TEG  yields 
the  same  power  output  but  with  a  higher  load  electrical  resistance. 
To  simplify  the  procedure  of  determining  mopt  we  have  developed  a 
performance  metric  based  on  the  rate  of  change  of  a  material’s  effi¬ 
ciency  as  a  function  of  the  material’s  working  temperature  gradi¬ 
ent  (/» =  ni  AT). 

As  depicted  in  Fig.  4,  Bi2Te3  with  a  fmax  of  4.57  x  10  2  K  1  per¬ 
forms  better  than  PbTe  which  reaches  a  pmax  of  1.86  x  10  2  K 
Note  that  pmax  defines  the  optimum  m  value  (mopt).  When  the  aver¬ 
age  a  value  is  considered  for  the  segmented  Bi2Te3-PbTe  (a)  TEG, 
pmax  reaches  a  peak  value  of  5.67  x  10  2  K  1  compared  to 
2.82  x  10  2  K  1  for  Bi2Te3-PbTe  (a).  In  terms  of  m,  the  optimum 
value  of  1.545  is  obtained  in  the  case  of  the  segmented  Bi2Te3- 
PbTe  (a)  with  the  TEG  efficiency  calculated  on  the  performance 
of  each  individual  material  segment.  As  expected,  the  maximum 
(Anax)  value  of  Bi2Te3-PbTe  (a)  of  2.82  x  10  2  K  1  resides  between 
the  flmax  values  of  Bi2Te3  and  PbTe  and  favorably  sustains  a  higher 
m  value  of  1.545.  The  aberrant  flmax  value  of  Bi2Te3-PbTe  (a)  is 
attributed  to  the  normalization  of  a. 


Electrical  Resistance  Ratio  (m=RL/RTEG) 

Fig.  4.  Ratio  of  load  (RJ  to  TEG  (Rtec)  electrical  resistance  as  a  function  of  the  TEG 
temperature-dependent  performance.  Solid  black  circles  denote  m„pt,  at  the 
maximum  /i  value,  for  each  material  or  design  configuration.  The  average  and  the 
collective  Seebeck  coefficient  cases  are  denoted  by  a  and  a,  respectively. 


4.2.  Integration  of  the  analytical  and  computational  models 


Due  to  the  interdependency  between  the  TEG  arms’  length  ( t ) 
and  the  Bi2Te3-PbTe  interface  temperature  (T/nter),  an  iterative  pro¬ 
cess  was  adopted  to  determine  the  latter  parameters.  Initially,  the 
boundary  condition  temperatures  TH,  Tc,  and  Tlnterwere  assumed  to 
be  623  K,  298  K,  and  475  K,  respectively.  Following  the  procedure 
outlined  in  Subsection  3.1  the  efficiency  and  geometry  of  the  TEG 
configurations  were  determined.  The  algorithm  for  obtaining  the 
boundary  condition  temperatures  TH,  Tc,  and  the  mean  Bi2Te3- 
PbTe  interface  temperature  (TJnter)  is  based  on  an  iterative  proce¬ 
dure  presented  in  Fig.  5. 


4.3.  TEG-heat  sink  thermal  resistance  matching 

As  a  prerequisite  for  achieving  the  maximum  TEG  power,  the 
TEG  thermal  resistance  and  the  heat  sink  and  heat  source  thermal 
resistance  should  match  [1,24],  To  investigate  this  assertion,  the 
performance  of  the  four  design  configurations  (Bi2Te3,  PbTe,  Bi2_ 
Te3-PbTe  (a),  Bi2Te3-PbTe  (a))  was  determined.  For  this  purpose 
the  TEG  thermal  resistance  (©teg)  to  the  heat  source  and  heat 
exchanger  thermal  resistance  (©Hx),  denoted  by  n  =  0tegI&hx, 
was  used.  Thermal  resistance  was  implemented  using  0teg,j  = 
%Gy/(KiEcw  -  ATf.c:j )  where  i  and  j  refer  to  the  n-  and  p-type  materi¬ 
als  of  Bi2Te3  and  PbTe.  The  latter  equation  was  used  to  determine 
the  thermal  power  under  no-load  (PH)  and  load  (PL)  conditions, 
for  a  range  of  n  values,  which  defines  the  efficiency  of  the  Bi2Te3- 
PbTe  system  (Eq.  (8)).  Contacts  between  the  TE  elements  were  as¬ 
sumed  to  be  ideal  therefore  not  negatively  affecting  the  thermal 
performance  of  the  TEG  configurations  examined  herein.  Normal¬ 
ized  changes  in  efficiency  (q)  were  used  to  identify  the  thermal 
resistance  matching  point  between  ©teg  and  0Hx  which  produced 
the  peak  TEG  performance,  that  is,  r\  =  1,  as  shown  in  Fig.  6. 

Examining  the  trend  of  matching  ©teg  and  ©hx,  as  depicted  in 
Fig.  6,  it  becomes  apparent  that  a  correlation  exists  between  ther¬ 
mal  resistance  and  the  temperature  gradient  (AT)  within  which 
each  design  configuration  operates.  Thermal  matching  (@jeg  =  - 
©hx)  for  Bi2Te3  occurs  at  5.06  KW  1  for  AT  =  153.2  K.  For  PbTe, 
when  AT  increases  to  171.4  K,  &teg=&hx  reduces  to  4KW  1 . 
Notably,  Bi2Te3-PbTe  (a)  and  Bi2Te3-PbTe  (a)  register  the 
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Fig.  5.  Iterative  process  of  determining  the  efficiency  and  geometry  of  segmented 
Bi2Te3-PbTe  TEGs  based  upon  the  computational  and  analytical  models. 

maximum  thermal  resistance  matching  point  (  &teg  =  &hx)  of 
12.62  KVT1  when  AT  increases  to  324.6  K. 

Compared  to  Bi2Te3  and  PbTe,  the  average  (a)  and  collective  (a) 
Bi2Te3-PbTe  designs  exhibit  a  considerably  higher  (60-68%)  source 
and  sink  thermal  resistance  matching  (&teg  =  &Hx)-  Apparently, 
matching  the  thermal  resistance  of  the  heat  source  and  heat  ex¬ 
changer  to  the  TEG  thermal  resistance  is  another  necessary  step  to¬ 
wards  optimizing  the  TEG  performance. 

4.4.  Segmented  and  average  seebeck  coefficient  Bi2Te3-PbTe  TECs 

A  simplified  way  to  derive  the  efficiency  of  a  segmented  TEG  is 
to  average  the  Seebeck  coefficient  (a)  values  of  the  constituent 


Thermal  Resistance  Matching  (0teg  =  0h*,  K-W  l) 


Fig.  6.  Normalized  changes  in  efficiency  as  a  function  of  the  thermal  resistance 
matching  between  the  TEG  ( G-nr.)  and  the  heat  source  and  heat  exchanger  ( 0Hx). 
Maximum  TEG  system  performance  is  attained  when  (~)Hx  matches  @rec,  i.e.,  *1  =  1. 


collective  (a)  segmented  TEGs.  The  reduced  material  efficiencies  for  a  given 
configuration  indicate  the  TEG  to  heat  source  and  heat  exchanger  thermal 

elements  over  the  TEG  operating  temperature  range.  Using  the 
average  Seebeck  coefficient  (Bi2Te3-PbTe),  Voc  can  be  obtained  for 
AT  =  324.6  K.  Subsequently,  m  (Bi2Te3-PbTe)  and  n  (Bi2Te3-PbTe) 
can  be  determined  from  the  methodology  outlined  in  Subsection 
3.1.  Eventually,  the  full  load  thermal  power  ( PH )  of  the  TEG,  upon 
which  the  TEG  efficiency  (ij)  is  based,  is  derived.  The  efficiency  of 
the  segmented  TEGs  is  obtained  by  considering  the  full  thermal 
power  produced  by  the  individual  thermoelements  of  Bi2Te3  and 
PbTe  over  their  respective  temperature  regime. 

Averaging  the  Seebeck  coefficient  (a)  values  of  Bi2Te3  and  PbTe 
within  the  entire  temperature  gradient  of  the  TEG,  of  324.6  K,  ap¬ 
pears  to  overestimate  the  Bi2Te3-PbTe  peak  efficiency  (rj.  Fig.  7). 
However,  when  the  efficiency  of  the  segmented  Bi2Te3-PbTe  con¬ 
figuration  is  derived  from  the  performance  of  the  individual  Bi2Te3 
and  PbTe  materials,  based  on  their  corresponding  a  values,  an »;  va¬ 
lue  of  5.29%  was  obtained  (Bi2Te3-PbTe  (a),  Fig.  7).  The  latter  value 
is  in  very  good  agreement  with  other  investigator’s  findings  [1,7], 
Matching  the  thermal  resistance  of  the  TEG  to  the  thermal  resis¬ 
tance  of  the  heat  source  and  the  heat  exchanger  (n  =  &tegI 
&hx  =  1).  as  depicted  in  Fig.  7,  results  in  an  increase  maximum  effi¬ 
ciency  of  74.3%  for  Bi2Te3  (a)  and  140.9%  for  Bi2Te3  (a)  relative  to 
n  =  3. 

Results  also  indicate  that  for  a  given  temperature  gradient  (AT) 
the  magnitude  of  the  Seebeck  coefficient  is  the  single  most  impor¬ 
tant  parameter  which  governs  efficiency.  In  other  words,  theoret¬ 
ical  performance  calculations  largely  reflect  both  the  theoretical 
methodology  and  the  accuracy  of  ct  values.  Clearly,  further 
progress  should  factor  in  the  temperature  dependent  variation 


4.5.  Geometry  of  the  TE  elements 

The  optimum  TEG  dimensions  which  maximize  the  TEG  device 
figure  of  merit  (and  efficiency)  are  derived  from  Eq.  (11).  Results 
for  Bi2Te3,  PbTe,  Bi2Te3-PbTe  (a),  and  Bi2Te3-PbTe  (a)  for  the  cor¬ 
responding  temperature  gradient  of  interest  and  Jmax,  yield  differ¬ 
ent  TE  lengths,  as  depicted  in  Fig.  8.  Of  particular  interest  is  the 
more  compact  size  of  the  segmented  Bi2Te3-PbTe  (a)  whose  length 
measures  0.152  cm  (or  4.6%  shorter)  in  comparison  to  Bi2Te3-PbTe 
(a)  which  has  a  length  of  0.159  cm.  Results  also  suggest  an 
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Fig.  8.  Length  of  the  thermoelement  (TE)  couples  of  the  Bi2Te3  and  PbTe  materials 
and  the  average  Bi2Te3-PbTe  (a)  and  collective  Bi2Te3-PbTe  (a)  configurations. 

inversely  proportional  correlation  between  the  TEG  arms’  lengths 
and  the  maximum  current  density. 

5.  Conclusions 

Confirming  experimental  measurements,  analytical  modeling 
results,  augmented  by  computational  simulations,  of  the  maxi¬ 
mum  efficiency  to  temperature  gradient  =  J|)  showed  that 
the  collective  Seebeck  coefficient  Bi2Te3-PbTe  (a)  design  can  sup¬ 
port  a  higher  electrical  load  in  relation  to  the  homogeneous  Bi2Te3 
and  PbTe  materials.  Averaging  the  Seebeck  coefficient  in  the  case  of 
Bi2Te3-PbTe  (a)  TEG  overestimates  pmax,  at  the  highest  m  value, 
compared  to  the  Bi2Te3,  PbTe,  and  Bi2Te3-PbTe  (£)  TEGs.  Optimum 
TEG  performance  (17)  in  terms  of  matching  the  TEG  to  the  heat 
source  and  heat  exchanger  thermal  resistance  ratio  revealed  that 
the  thermal  resistance  of  the  segmented  TEG  system  increases 
with  the  temperature  gradient  of  the  TEG.  Matching  ©teg  to  &Hx 
is  a  prerequisite  for  maximizing  the  TEG  system  efficiency  and  con¬ 
trolling  the  size  of  the  heat  exchanger  which,  in  turn,  governs  the 
physical  dimensions  of  the  TEG.  Actively  cooled  heat  exchangers, 
as  our  results  demonstrate,  constitute  an  attractive  option  for  high 
thermal  energy  density  dissipation  and  compact  TEGs. 

Computational  results  provided  exquisite  details  of  the  temper¬ 
ature  distribution  across  the  TEG  legs  which,  in  turn,  helped  raise 
the  accuracy  of  the  analytical  results.  In  closing,  the  hybrid  compu¬ 
tational-analytical  modeling  approach  achieves  a  very  good  level 
of  accuracy  as  part  of  a  time  efficient  framework. 
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